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Search for Contact Interactions in Dilepton Events from pp Collisions 
at y/s = 7 TeV with the ATLAS Detector 

The ATLAS Collaboration 

This Letter presents a search for contact interactions in the dielectron and dimuon channels using 
data from proton-proton collisions produced by the LHC at y/s = 7 TeV and recorded by the 
ATLAS detector. The data sample, collected in 2011, corresponds to an integrated luminosity of 
1.08 and 1.21 fb" 1 in the e+e" and H + fj, channels, respectively. No significant deviations from the 
standard model are observed. Using a Bayesian approach with a prior flat in 1/A 2 , the following 
95% CL lower limits are placed on the energy scale of Uqq contact interactions: A~ > 10.1 TeV 
(A+ > 9.4 TeV) in the electron channel and A" > 8.0 TeV (A+ > 7.0 TeV) in the muon channel 
for constructive (destructive) interference in the left-left isoscalar contact interaction model. Limits 
are also provided for a prior flat in 1/A 4 . 



INTRODUCTION 

A wide range of new physics phenomena can produce 
modifications to the dilepton mass spectra predicted by 
the standard model (SM) such as quark/lepton compos- 
iteness, extra dimensions, and new gauge bosons. The 
predicted form of these deviations is often either a reso- 
nance or an excess in the number of events in the spectra 
at high mass. This Letter reports on a search for such 
an excess in dilepton events produced in proton-proton 
collisions at the LHC [l[ . An interpretation of these data 
in the context of contact interactions (CI) is presented, 
including the first limits with the ATLAS detector in the 
dielectron channel and an update of the search performed 
using 2010 data in the dimuon channel [2|. A separate pa- 
per describes the search for new heavy resonances in the 
dilepton mass spectra performed using the same ATLAS 
dataset 0. 

If quarks and leptons are composite, with at least 
one common constituent, the interaction of these con- 
stituents would likely be manifested through an effec- 
tive four-fermion contact interaction at energies well be- 
low the compositeness scale. Such a contact interaction 
could also describe a new interaction with a messenger 
too heavy for direct observation at the LHC, in analogy 
with Fermi's nuclear /? decay theory [4j. 

The Lagrangian for a general contact interaction has 
the form [5| 

+ VRR jPrIh^R^PrI^R C 1 ) 
+2r] LR iPlJ^l iPr^iPr } , 

where g is a coupling constant chosen to obey g 2 /47r = 1, 
A is the contact interaction scale, which in the context 
of compositeness models is the energy scale below which 
fcrmion constituents are bound, and ipL,R are left-handed 
and right-handed fermion fields, respectively. The pa- 
rameters rjij , where i and j are L or R, define the chiral 
structure (left or right) of the new interaction. Specific 



models are constructed by setting different combinations 
of these parameters to assume values of —1, or +1. 
The addition of this contact interaction term to the SM 
Lagrangian alters the Drell-Yan (DY) production cross 
section (qq — > Z/j* — > The largest deviations, 

either constructive or destructive, are expected at high 
dilepton invariant mass and are determined by the scale 
A and the sign of the parameter rjij. This analysis in- 
terprets the data in the context of the left-left isoscalar 
model (LLIM), which is commonly used as a benchmark 
for contact interaction searches [6] . The LLIM is defined 
by setting rj LL = ±1 and rj RR = r} L R = 0. 

With the introduction of a contact interaction, the dif- 
ferential cross section for the process qq —> £ + £~ can be 
written 

da _ dapY ^^ Fi(m u ) | F c (m u ) , 
dmu dmu A 2 A 4 

where mu is the final-state dilepton mass. The expres- 
sion above includes a SM DY term, as well as DY-CI in- 
terference (Fi) and pure contact interaction (Fq) terms 
(see Ref . [7| for the full form of this expression) . At the 
largest A values to which this analysis is sensitive, both 
interference and pure contact interaction terms play a 
significant role. For example, at dilepton masses greater 
than 300 GeV and A = 9 TeV, the magnitude of the inter- 
ference term is about 1.5 times that of the pure contact 
interaction term. 

The present analysis focuses on identifying a broad de- 
viation from the SM dilepton mass spectra, which are 
expected to be dominated by the DY process. Current 
experimental bounds on A (see below) indicate any devi- 
ation from a new interaction would appear at masses well 
above the Z boson peak. Consequently, the search region 
is restricted to dilepton masses above 150 GeV. The anal- 
ysis exploits the high pp collision energy of the LHC and 
the capabilities of the ATLAS detector to identify and 
reconstruct electrons and muons at high momentum. 

Previous searches for contact interactions have been 
carried out in neutrino scattering 18], as well as at 
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electron-positron |9h13I1 . electron-proton [lj, [l5|, and 
hadron colliders [16l424| . In the case of eeqq contact 
interactions, the best limits in the LLIM for all quark 
flavors come from e + e~ experiments with A~ > 7.2 TeV 
and A+ > 12.9 TeV @ at 95% confidence level (CL) for 
= —1 and +1, respectively. These limits assume 
that contact interactions of electrons with all quark fla- 
vors are of the same strength. Best limits set in the 
specific case of first generation quarks are A~ > 9.1 TeV 
and A+ > 8.6 TeV |l3j at 95% CL. In the case of fifiqq 
contact interactions, the best limits are A~ > 4.9 TeV 
and A+ > 4.5 TeV from the ATLAS analysis of the 2010 
data 0). 



ATLAS DETECTOR AND DATA SAMPLE 

ATLAS is a multipurpose particle detector [25| . It con- 
sists of an inner tracking detector surrounded by a 2 T 
superconducting solenoid, electromagnetic and hadronic 
calorimeters, and a muon spectrometer with a toroidal 
magnetic field. Charged particle tracks are reconstructed 
using the inner detector, which comprises a silicon pixel 
detector, a silicon-strip tracker, and a transition radi- 
ation tracker, covering the pseudorapidity range \r]\ < 
2.5 [2(|. A hermetic calorimeter, which covers \r)\ < 
4.9, surrounds the superconducting solenoid. The liquid- 
argon electromagnetic calorimeter, which plays an im- 
portant role in electron identification and measurement, 
is finely segmented, with readout granularity (77, <f) vary- 
ing by layer and cells as small as 0.025 x 0.025 ex- 
tending to 1 77 1 < 2.5, to provide excellent energy and 
position resolution. The electron energy resolution is 
dominated at high energy by a constant term equal to 
1.2% in the barrel (\rj\ < 1.37) and 1.8% in the endcaps 
(1.52 < |7y| < 2.47). Hadron calorimetry is provided by 
an iron-scintillator tile calorimeter in the central rapid- 
ity range \rj\ < 1.7 and a liquid- argon calorimeter in the 
rapidity range 1.5 < I77I < 4.9. Another key detector 
component for this analysis is the muon spectrometer, 
which is designed to identify muons and measure their 
momenta with high accuracy. The currently achieved 
resolution for momenta transverse to the beam line (pr) 
of 1 TeV ranges from 15% (central) to 44% (for \r)\ > 2). 
The muon system comprises three toroidal magnet sys- 
tems, a trigger system consisting of resistive plate cham- 
bers in the barrel and thin-gap chambers in the endcaps, 
providing triggering capability up to \r]\ = 2.4, and a 
set of precision monitored drift tubes and cathode strip 
chambers in the region \r)\ < 2.7. 

The data sample for this analysis was collected during 
LHC operation in the first half of 2011 and corresponds to 
a total integrated luminosity of 1.08 and 1.21 fb _1 in the 
e + e~ and fi >fj,~ channels, respectively. It was collected 
with stable beam conditions and an operational inner de- 
tector. For the electron (muon) channel, the calorimeter 



(muon spectrometer) was also required to be operational. 
Events were selected by requiring that they pass the sin- 
gle electron (muon) trigger with a transverse momentum 
Pt threshold of 20 (22) GeV. This analysis follows the 
same event selection as the search for new heavy reso- 
nances. A summary is provided below, a more complete 
description can be found in Ref. 



SIGNAL AND BACKGROUND MODELING 

This analysis looks for deviations from the expected 
SM dilepton spectra. The largest SM contribution comes 
from DY followed by semileptonic decay of tt pairs, 
electroweak diboson production (WW, WZ, and ZZ), 
and production of jets in association with a W boson 
(W+jets). In addition, multi-jet production (QCD) is 
a significant background in the electron channel. With 
the exception of QCD, Monte Carlo (MC) simulation was 
used to model these backgrounds. 

DY events were generated with Pythia 6.421 [53] 
and Mrst2007 LO* parton distribution functions 
(PDFs) [H. Signal DY+CI samples in the LLIM were 
generated with the same version of Pythia for the full 
dilepton differential cross section as shown in Eq. ([2]l. 
This ensured that the interference term Fj was properly 
included. All quark flavors contributed to the contact 
interaction in these signal samples. Diboson processes 
were produced with Herwig 6.510 j2§| using Mrst2007 
LO* PDFs. The VF+jets background was generated with 
Alpgen [13] and Cteq6l1 J3p PDFs, and the tt back- 
ground with Mc@NLO 3.41 [33 and Cteq6.6 (33j PDFs. 
For the latter two, Jimmy 4.31 [34j was used to describe 
multiple parton interactions and Herwig to describe the 
remaining underlying event and parton showers. Pho- 
tos [sH was used to handle the final-state photon radia- 
tion for all MC samples. Furthermore, higher order QCD 
corrections were implemented via a mass-dependent K- 
factor defined as the ratio between the next-to-next-to- 
leading order (NNLO) Z/j* cross section, calculated us- 
ing Phozpr [S] and Mstw2008 PDFs and the LO 
cross section. This QCD A'-factor was applied to both 
DY and DY+CI samples. Likewise, DY and DY+CI 
samples are corrected with a mass-dependent if-factor 
accounting for higher-order electroweak corrections aris- 
ing from virtual heavy gauge boson loops that are calcu- 
lated using Horace [38(. Finally, the generated samples 
were processed through a full simulation of the ATLAS 
detector based on the Geant 4 package [40j | . 

For both channels, the QCD multi-jet background 
is evaluated from data due to poor modeling and low 
MC statistics. In the electron channel, a reversed elec- 
tron identification technique is used to select a sample 
of events in which both electrons fail a subset of the 
electron identification criteria (see further discussion be- 
low). This sample is then used to determine the shape 
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of the QCD background as a function of dielectron in- 
variant mass. This template shape and the sum of the 
DY, dibosons, tt, and W+jets backgrounds normalized 
by their cross sections (including higher order correc- 
tions) are fitted to the observed dielectron mass distri- 
bution in the range between 70 and 200 GeV to deter- 
mine the normalization of the QCD contribution. The 
above QCD background estimate is cross-checked with 
two other methods described in Ref. [H in order to de- 
termine its systematic uncertainty. In particular, these 
cross-checks set bounds on the potential bias in the QCD 
mass spectrum introduced by the reversed identification 
technique. In the muon channel, the QCD background 
is much smaller and is also evaluated from data. A re- 
verse isolation method is utilized: a QCD sample is se- 
lected from data by requiring two non-isolated muons 
with 0.1 < Ep T (i? < 0.3)/p T (» < 1.0, where Sp T 
(R < 0.3) is the sum of the pt of the tracks in a cone 
of R = yj A77 2 + Acf) 2 < 0.3 around the direction of the 
muon. The normalization for this sample is obtained 
from the ratio of isolated to non-isolated dimuon events 
in QCD cc and bb MC, where the isolation requirement 
is YiPt(R < 0.3)/pt(m) < 0.05. Muons from light hadron 
decays are not a significant source of background at the 
high momenta relevant to this analysis. 



EVENT SELECTION 

Events passing the trigger selection described above 
are required to have a pair of either electrons or muons 
with pt greater than 25 GeV to ensure maximal trigger 
efficiency. To reject cosmic ray events and beam halo 
background, events are required to have a reconstructed 
vertex with at least three charged particle tracks with 
Pt > 0.4 GeV. If several such vertices are found, the ver- 
tex with the largest ^p T is selected as the primary vertex 
of the event, where the sum is over all charged particles 
associated with the given vertex. Electron candidates are 
confined in \t)\ < 2.47, with the detector crack region 1.37 
< 1 77 1 < 1.52 excluded because of degraded energy reso- 
lution. Muon candidates are required to be within the 
inner detector acceptance. 

Electron candidates are formed from clusters of cells 
reconstructed in the electromagnetic calorimeter. Identi- 
fication criteria on the transverse shower shape, the lon- 
gitudinal leakage into the hadronic calorimeter, and the 
association to an inner detector track are applied to the 
cluster to satisfy the medium electron definition 4l( . The 
electron energy is obtained from the calorimeter mea- 
surements and its direction from the associated inner 
detector track. A hit in the first layer of the pixel de- 
tector is required (if an active pixel layer is traversed) to 
suppress background from photon conversions. Further 
QCD jet background suppression is achieved by demand- 
ing the highest pt electron in the event to be isolated. 



To that effect, the sum of the calorimeter transverse mo- 
menta around the electron direction Y<pt(R < 0.2) must 
be less than 7 GeV. The core of the electron energy depo- 
sition is excluded and the sum is corrected for transverse 
shower leakage and pile-up from additional pp collisions. 
In addition, the two electron candidates are not required 
to have opposite charge because of possible charge mis- 
identification either due to bremsstrahlung or to the lim- 
ited momentum resolution of the inner detector at very 
high pt- If the event contains more than two selected 
electrons, the two electrons with the highest pt are cho- 
sen. For these selection criteria, the overall event accep- 
tance for signal events has a small dependence on the 
dielectron mass above 500 GeV and a value of approxi- 
mately 65% at 1 TeV. 

Muon candidates are required to be of opposite charge 
and are reconstructed independently in both the inner 
detector and the muon spectrometer. The momentum 
is taken from a combined fit to the measurements from 
both subsystems. To obtain optimal momentum resolu- 
tion and accurate modeling by the simulation, the muon 
candidates are required to satisfy the following require- 
ments in the muon spectrometer: have at least three hits 
in each of the inner, middle, and outer detectors, and 
at least one hit in the non-bending xy plane. To sup- 
press background from cosmic rays, requirements are im- 
posed on the muon impact parameter and primary vertex 
(PV): transverse impact parameter |do| < 0.2 mm, z co- 
ordinate with respect to the PV \zq — Zpvi < 1 mm, 
and z position of the PV |zpy| < 200 mm. Muons 
are also required to be isolated to reduce background 
from jets: T,pt(R < 0.3)/pt(m) < 0.05, as explained 
in the previous section. If more than one opposite-sign 
muon pair is found in an event, the pair with the largest 
Pt(a* + ) + Pt(m~) is chosen. The overall event accep- 
tance for signal events has only a weak dependence on 
the dimuon mass with a value of approximately 40% at 
1 TeV. Stringent requirements on the presence of hits in 
all three layers of the muon spectrometer and the limited 
three-layer geometrical coverage are the primary reason 
for the lower acceptance relative to the electron channel. 

Extensive comparisons between data and MC simula- 
tion were performed at the level of single-lepton distri- 
butions to confirm that the simulation reproduces the 
selected data well, especially at high momentum. 

Figure [1] displays the dielectron and dimuon mass spec- 
tra for all selected events with invariant mass greater 
than 70 GeV. The expected event yields for the differ- 
ent processes are obtained by first normalizing each MC 
process by its cross section (including higher order cor- 
rections) and then normalizing the total MC event yield 
plus the data-derived QCD background to the data in the 
Z peak region (dilepton mass between 70 and 110 GeV). 
Good agreement is observed between the data and the 
SM prediction over the whole dilepton mass range. A 
quantitative comparison is provided in Tables HI and [TT1 A 
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FIG. 1. Dielectron (top) and dimuon (bottom) invariant mass 
distributions for data (points) and Monte Carlo simulation 
(histograms). The open histograms correspond to the distri- 
butions expected in the presence of contact interactions with 
different values of A for both constructive (solid histograms) 
and destructive (dashed histograms) interference. 

slight excess of events observed at high dimuon mass is 
consistent with a statistical fluctuation. The most signif- 
icant deviation in the number of dimuon events occurs for 
events with mass greater than 800 GeV. In this region, 
the Poisson probability for observing 5 or more events 
where 2.1 are expected is 6.2%. The muon tracks in the 
five data events were inspected in detail and no problem 
was found. 



SYSTEMATIC UNCERTAINTIES 

Since the MC event yields are normalized to the num- 
ber of events observed in the Z peak region, only mass- 
dependent systematic uncertainties need to be consid- 
ered, except for a 5% overall uncertainty in the knowl- 
edge of the Z/'y* cross section in the normalization re- 
gion. This overall uncertainty is required since the cross 
section change due to the new physics is defined with re- 



spect to the SM cross section. The dominant uncertain- 
ties are of theoretical origin but experimental sources are 
also considered. 

Theoretical uncertainties in the predicted event yields 
arise from the limited knowledge of PDFs, as, QCD and 
electroweak K- factors, and Z/j* cross section. The finite 
available MC statistics are also taken into account. The 
uncertainty due to the PDF and as is estimated using the 
MSTW2008 PDF eigenvector set and additional PDFs 
corresponding to variations in as- The resulting effect 
is about 4% at the Z pole growing with increasing mass 
to 10% at 1.5 TeV. Uncertainties due to QCD and elec- 
troweak if-factors are estimated to grow from 0.3% and 
0.4% at the Z pole to 3% and 4.5% at 1.5 TeV, respec- 
tively, for both electron and muon channels. Estimates 
for the QCD X-factor are obtained by varying the renor- 
malization and factorization scales independently by fac- 
tors of two, then adding the impact of those variations 
linearly. The uncertainty in the electroweak isT-factor is 
evaluated from the effect of neglecting real boson emis- 
sion, varying the electroweak scheme definitions as im- 
plemented in Pythia and Horace, as well as of the ef- 
fect of higher order electroweak and O(aas) corrections. 
For the electron channel, the QCD background estimate 
is subject to an uncertainty derived from a comparison 
with different background estimate methods (see discus- 
sion above). For the muon channel, the QCD background 
uncertainty is negligible. 

Experimental uncertainties originate from the en- 
ergy/momentum resolution, as well as the trigger, re- 
construction and identification efficiencies. In the elec- 
tron channel, the uncertainty in the constant term, which 
dominates the energy resolution at high energy, has a 
negligible impact on the analysis. Knowledge of the en- 
ergy scale also has a negligible effect. The electron re- 
construction and identification uncertainty results in a 
1.5% effect, which is estimated by studying the impact 
of the isolation requirement on the dielectron mass distri- 
bution. In the muon channel, the momentum resolution 
is dominated by the quality of the muon spectrometer 
alignment. The uncertainty in the alignment is evalu- 
ated directly from dedicated toroid field-off runs and re- 
dundant momentum measurements in overlapping small 
and large chambers. These experimental uncertainties 
are found to have minimal impact on the dimuon mass 
distribution. Finally, a systematic error growing from 
0.3% at the Z pole to 4.5% at 1.5 TeV is assigned to 
the muon reconstruction efficiency to account conserva- 
tively for its small px dependence due to occasional large 
energy loss from bremsstrahlung. 

STATISTICAL ANALYSIS 

The data analysis proceeds with a Bayesian method 
to compare the observed event yields with the expected 
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TABLE I. Expected and observed numbers of events in the electron channel. The total expected yield is normalized to the 
data in the Z peak control region between 70 and 110 GeV. The errors quoted originate from both systematic uncertainties 
and limited MC statistics, except the error on the total expected yield in the normalization region which is given by the square 
root of the number of observed events. 



m cc [GeV] 


70-110 


110-130 


130-150 


150-170 


170-200 


200-240 


240-300 


DY 


258482 ± 410 


3185 ± 110 


1183 ± 46 


608 ± 28 


473 ± 24 


312 ± 18 


196 ± 9 


ti 


218 ± 36 


87 ± 7 


64 ± 5 


51 ± 2 


51 ± 2 


37 ± 2 


30 ± 1 


Dibosons 


368 ± 19 


31 ± 2 


24 ± 2 


15 ± 1 


16 ± 1 


14 ± 1 


8 ± 1 


W+jcts 


150 ± 100 


57 ± 11 


40 ± 9 


27 ± 6 


26 ± 6 


20 ± 5 


14 ± 4 


QCD 


332 ± 60 


80 ± 42 


50 ± 20 


32 ± 7 


29 ± 7 


19 ± 15 


12 ± 9 


Total 


259550 ± 510 


3440 ± 120 


1361 ± 50 


733 ± 30 


595 ± 25 


401 ± 24 


260 ± 13 


Data 


259550 


3419 


1362 


758 


578 


405 


256 




m ee [GeV] 


300-400 


400-550 


550-800 


800-1200 


1200-1800 


1800-3000 




DY 


105.0 ± 5.0 


41.0 ± 2.2 


12.8 ± 0.8 


2.5 ± 0.2 


0.29 ± 0.05 


< 0.05 




ti 


14.9 ± 0.8 


4.5 ± 0.2 


1.0 ± 0.1 


0.10 ± 0.02 


< 0.05 


< 0.05 




Dibosons 


7.5 ± 1.1 


2.1 ± 0.4 


1.0 ± 0.3 


0.3 ± 0.1 


< 0.05 


< 0.05 




W+]ets 


9.0 ± 3.2 


3.5 ± 1.6 


1.0 ± 0.7 


0.2 ± 0.3 


< 0.05 


< 0.05 




QCD 


5.5 ± 4.4 


1.5 ± 1.2 


0.3 ± 0.2 


< 0.05 


< 0.05 


< 0.05 




Total 


141.9 ± 8.0 


52.6 ± 3.0 


16.1 ± 1.1 


3.0 ± 0.4 


0.33 ± 0.05 


< 0.05 




Data 


147 


48 


17 


3 











TABLE II. Expected and observed numbers of events in the dimuon channel. The total expected yield is normalized to the 
data in the Z peak control region between 70 and 110 GeV. The errors quoted originate from both systematic uncertainties 
and limited MC statistics, except the error on the total expected yield in the normalization region which is given by the square 
root of the number of observed events. 



winn [GeV] 


70-110 


110-130 


130-150 


150-170 


170-200 


200-240 


240-300 


DY 


236405 ± 320 


3133 ± 90 


1076 ± 36 


548 ± 22 


417 ± 18 


249 ± 13 


153 ± 7 


ti 


193 ± 21 


70 ± 9 


51 ± 7 


34 ± 4 


38 ± 4 


30 ± 3 


21 ± 2 


Diboson 


307 ± 16 


25 ± 2 


19 ± 2 


13 ± 2 


12 ± 1 


10 ± 1 


8 ± 1 


W+jets 


1 ± 1 


< 0.5 


< 0.5 


< 0.5 


< 0.5 


< 0.5 


< 0.5 


QCD 


1 ± 1 


< 0.5 


< 0.5 


< 0.5 


< 0.5 


< 0.5 


< 0.5 


Total 


236908 ± 490 


3229 ± 90 


1147 ± 37 


595 ± 22 


467 ± 19 


290 ± 13 


182 ± 8 


Data 


236908 


3211 


1132 


621 


443 


279 


195 






300-400 


400-550 


550-800 


800-1200 


1200-1800 


1800-3000 




DY 


80.8 ± 3.9 


31.0 ± 1.7 


9.2 ± 0.6 


1.8 ± 0.2 


0.22 ± 0.04 


< 0.05 




ti 


11.7 ± 1.2 


3.5 ± 0.3 


0.7 ± 0.1 


0.06 ± 0.02 


< 0.05 


< 0.05 




Diboson 


6.7 ± 1.1 


1.0 ± 0.4 


0.7 ± 0.3 


< 0.05 


< 0.05 


< 0.05 




W+jots 


< 0.05 


< 0.05 


< 0.05 


< 0.05 


< 0.05 


< 0.05 




QCD 


< 0.05 


< 0.05 


< 0.05 


< 0.05 


< 0.05 


< 0.05 




Total 


99.3 ± 4.2 


35.5 ± 1.8 


10.6 ± 0.7 


1.9 ± 0.2 


0.22 ± 0.04 


< 0.05 




Data 


83 


39 


12 
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yields for a range of different contact interaction model 
parameters. Specifically, the number \x of expected 
events in each of the mass bins defined in Tables U and UU 
is 

/i = n DY +ci(#, v) + n non _ DY bg(^) , (3) 

where nDY+ci(0> v) is the number of events predicted by 
the Pythia DY+CI MC for a particular choice of con- 
tact interaction model parameter 6, ri non _DY bg(^) is the 
number of non-DY background events, and v represents 
the set of Gaussian nuisance parameters that account for 
systematic uncertainties in these numbers as discussed 



above. The parameter 9 corresponds to a choice of energy 
scale A and interference parameter tjll. The complete 
set of /x values used in this analysis is shown in Tables IIIII 
and IIVI for the electron and muon channels, respectively. 
For each mass bin, a second order polynomial is used to 
model the dependence of /x on 1/A 2 . 

The likelihood of observing a set of n events in TV in- 
variant mass bins is given by a product of Poisson prob- 
abilities for each mass bin k: 

£(n | M) = J] (4) 

k—1 



6 



TABLE III. Expected numbers of events in the signal region of the analysis for various contact interaction scales with con- 
structive (A - ) and destructive (A + ) interference in the electron channel. The errors quoted originate from both systematic 
uncertainties and limited MC statistics. 





:e [GcV] 


150-170 


170-200 


200-240 


240-300 


300-400 


A - 


= 3 TeV 


785 ± 29 


649 ± 26 


467 ± 22 


383 ± 19 


343 ± 12 


A~ 


= 4 TeV 


781 ± 28 


647 ± 26 


437 ± 21 


326 ± 17 


223 ± 7 


A~ 


= 5 TeV 


734 ± 27 


612 ± 24 


405 ± 19 


298 ± 16 


181 ± 6 


A~ 


= 7 TeV 


691 ± 26 


638 ± 25 


406 ± 19 


259 ± 15 


163 ± 5 


A" 


= 12 TeV 


721 ± 26 


604 ± 24 


336 ± 17 


234 ± 14 


149 ± 5 


A+ 


= 3 TeV 


770 ± 28 


642 ± 24 


424 ± 20 


331 ± 17 


269 ± 9 


A+ 


= 4 TeV 


745 ± 27 


591 ± 23 


385 ± 19 


277 ± 16 


166 ± 5 


A+ 


= 5 TeV 


702 ± 25 


607 ± 23 


350 ± 17 


258 ± 15 


151 ± 5 


A+ 


= 7 TeV 


672 ± 25 


600 ± 23 


399 ± 19 


251 ± 14 


142 ± 5 


A+ 


= 12 TeV 


749 ± 27 


593 ± 23 


403 ± 19 


274 ± 15 


137.9 ± 4.4 





;e [GcV] 


400-550 


550-800 


800-1200 


1200-1800 


1800-3000 


A~ 


= 3 TeV 


286 ± 11 


269 ± 12 


207 ± 11 


112 ± 8 


30.3 ± 3.4 


A~ 


= 4 TeV 


132 ± 5 


109.1 ± 4.9 


80.5 ± 4.1 


29.8 ± 2.3 


10.9 ± 1.3 


A~ 


= 5 TeV 


82.7 ± 3.6 


57.3 ± 2.6 


35.5 ± 1.9 


14.8 ± 1.1 


3.9 ± 0.5 


A" 


= 7 TeV 


68.3 ± 2.9 


29.0 ± 1.4 


11.2 ± 0.7 


4.0 ± 0.4 


1.08 ± 0.17 


Ar 


= 12 TeV 


57.3 ± 2.6 


18.8 ± 1.1 


5.0 ± 0.4 


1.00 ± 0.13 


0.15 ± 0.05 


A+ 


= 3 TeV 


215 ± 8 


239 ± 11 


185 ± 10 


107 ± 7 


28.4 ± 3.2 


A+ 


= 4 TeV 


100.9 ± 3.8 


78.1 ± 3.6 


60.7 ± 3.2 


31.5 ± 2.2 


8.7 ± 1.0 


A+ 


= 5 TeV 


64.4 ± 2.7 


36.2 ± 1.8 


25.4 ± 1.3 


12.7 ± 0.9 


3.5 ± 0.4 


A+ 


= 7 TeV 


56.3 ± 2.5 


20.3 ± 1.1 


7.7 ± 0.5 


3.58 ± 0.28 


0.83 ± 0.12 


A+ 


= 12 TeV 


52.4 ± 2.4 


14.4 ± 0.9 


3.29 ± 0.28 


0.46 ± 0.08 


0.075 ± 0.026 



According to Bayes' theorem, the posterior probability 
for the parameter 6 given n observed events is 

V(6 | ft) = ±£ M (n | 6)P{0), (5) 

where Z is a normalization constant and the marginal- 
ized likelihood Cm corresponds to the likelihood after 
all nuisance parameters have been integrated out. This 
integration is performed assuming that the nuisance pa- 
rameters are correlated across all mass bins and that they 
affect both signal and background expectations, except 
for the electroweak iif-factor that only affects the DY and 
DY+CI components. The prior probability P(8) is cho- 
sen to be flat in 1 / A 2 , motivated by the form of Eq. © . 
The 95% CL limit is then obtained by finding the value 
Ai im satisfying f£ Um V{6 | ft) dJd = 0.95, where 9 = I/A 2 . 
The above calculations have been performed with the 
Bayesian Analysis Toolkit (BAT) [IH . 



RESULTS 

To test the consistency between the data and the 
standard model, a likelihood ratio test was performed 
by producing a set of 1000 SM-like pseudoexperiments 
and comparing the likelihood ratio between the sig- 
nal+background and pure background hypotheses ob- 
tained in the data to the results of the pseudoexperi- 
ments. The signal+background likelihood is evaluated at 



the A value that maximizes it. The derived p-value, cor- 
responding to the probability of observing a fluctuation 
in the pseudoexperiments that is at least as signal-like as 
that seen in the data (i.e. with a maximum likelihood 
ratio greater or equal to that obtained in the data) , is es- 
timated to be 39% (79%) in the electron channel and 21% 
(5%) in the muon channel for constructive (destructive) 
interference. These values indicate that there is no sig- 
nificant evidence for contact interactions in the analyzed 
data and thus limits are set on the contact interaction 
scale A. 

Using the Bayesian method described above, the ex- 
pected 95% CL lower limit values on the energy scale A 
are found to be 9.6 ± 1.0 TeV for constructive interference 
(tjll = — 1) an d 9.3± 1.0 TeV for destructive interference 
(t]ll = +1) in the electron channel. The corresponding 
expected limits in the muon channel are 8.9±0.9 TeV and 
8.6±0.9 TeV. The quoted uncertainties correspond to the 
68% range of limits surrounding the median value of all 
limits obtained with a set of 1000 pseudoexperiments. 
Stronger limits are expected in the electron channel than 
in the muon channel due to the significantly larger ac- 
ceptance for the dielectron selection. 

The observed limits (at 95% CL) are A" > 10.1 TeV 
(A+ > 9.4 TeV) in the electron channel and A" > 8.0 TeV 
(A + > 7.0 TeV) in the muon channel for constructive 
(destructive) interference. These limits are summarized 
in Table El 

If instead of choosing the prior to be flat in 1/A 2 , it 
is selected to be flat in 1/A 4 to match the form of the 
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TABLE IV. Expected numbers of events in the signal region of the analysis for various contact interaction scales with con- 
structive (A - ) and destructive (A + ) interference in the muon channel. The errors quoted originate from both systematic 
uncertainties and limited MC statistics. 





[GcV] 


150-170 


170-200 


200-240 


240-300 


300-400 


A" = 


3 TcV 


638 ± 28 


547 ± 26 


371 ± 22 


285 ± 


19 


263 ± 13 


A" = 


4 TcV 


618 ± 27 


513 ± 24 


287 ± 18 


228 ± 


15 


163 ± 7 


A" = 


5 TcV 


572 ± 26 


478 ± 23 


357 ± 20 


206 ± 


14 


131 ± 6 


A" = 


7 TcV 


571 ± 25 


496 ± 23 


294 ± 18 


187 ± 


14 


113 ± 5 


A" = 


12 TcV 


606 ± 26 


441 ± 22 


252 ± 16 


191 ± 


14 


100.0 ± 3.8 


A+ = 


3 TcV 


602 ± 26 


417 ± 21 


332 ± 19 


205 ± 


15 


186 ± 10 


A+ = 


4 TcV 


575 ± 25 


456 ± 22 


286 ± 17 


182 ± 


13 


112 ± 5 


A+ = 


5 TcV 


554 ± 25 


483 ± 23 


289 ± 17 


167 ± 


12 


102.0 ± 4.0 


A+ = 


7 TcV 


557 ± 24 


435 ± 21 


292 ± 18 


196 ± 


14 


102.0 ± 4.4 


A+ = 


12 TeV 


576 ± 25 


421 ± 21 


256 ± 16 


186 ± 


13 


100.0 ± 3.9 





[GcV] 


400-550 


550-800 


800-1200 


1200-1800 


1800-3000 


A" = 


3 TcV 


192 ± 11 


185 ± 12 


151 ± 12 


60 ± 9 


18 ± 7 


A" 


4 TcV 


97 ± 5 


73 ± 5 


50.3 ± 4.0 


17.9 ± 2.9 


6.8 ± 2.4 


\ 


5 TcV 


62.8 ± 3.3 


37.5 ± 2.2 


21.8 ± 1.8 


8.7 ± 1.3 


2.8 ± 1.0 


\ = 


7 TcV 


45.9 ± 2.4 


19.5 ± 1.2 


7.9 ± 0.7 


3.0 ± 0.5 


0.59 ± 0.38 


A" = 


12 TcV 


38.0 ± 2.0 


14.0 ± 0.9 


2.90 ± 0.26 


0.91 ± 0.24 


0.11 ± 0.14 


A+ = 


3 TcV 


152 ± 9 


153 ± 11 


131 ± 11 


63 ± 8 


19 ± 6 


A+ = 


4 TcV 


68.2 ± 3.7 


54.7 ± 3.7 


42.8 ± 3.4 


22.2 ± 2.7 


7.4 ± 2.1 


A+ = 


5 TcV 


51.0 ± 2.5 


24.9 ± 1.6 


15.7 ± 1.3 


8.0 ± 1.0 


2.9 ± 0.8 


A+ = 


7 TcV 


33.9 ± 1.8 


13.3 ± 0.8 


5.04 ± 0.43 


1.85 ± 0.32 


0.63 ± 0.30 


A+ = 


12 TcV 


34.0 ± 1.8 


10.1 ± 0.6 


1.73 ± 0.18 


0.25 ± 0.12 


0.07 ± 0.10 



pure CI term in Eq. the observed limit in the elec- 
tron channel becomes weaker by 0.9 TeV (0.8 TeV) for 
constructive (destructive) interference. The correspond- 
ing shift to lower values is 0.4 TeV (0.3 TeV) in the muon 
channel, see Table [V] 

The quoted limits have been obtained with electroweak 
corrections applied to both DY and DY+CI samples for 
consistency, although part of the electroweak corrections 
cannot be computed reliably due to the unknown new 
physics represented by the contact interaction. However, 
this particular choice does not affect the observed limits 
significantly and results in slightly more conservative lim- 
its. The limits obtained without electroweak corrections 
improve by less than 0.1 TeV. 

Finally, a limit is set on the combination of the elec- 
tron and muon channels, assuming lepton universality, 
by computing a combined posterior probability for the 
two channels. The following sources of systematic un- 
certainty are treated as fully correlated between the two 
channels: PDF and as, QCD and electroweak if- factors, 
and Z/"f* cross section for normalization. The resulting 
combined limits are A" > 10.2 TeV and A+ > 8.8 TeV 
for the 1/A 2 prior. Table \V\ summarizes all limits for the 
two priors considered in this study and shows that the 
combined limit is weaker than the electron channel limit 
in the case of destructive interference. This is due to the 
slight excess of events observed in the muon channel for 
masses greater than 800 GeV, an excess that gives an 
increase in the likelihood that is stronger for destructive 
interference than it is for constructive interference. 



TABLE V. Expected and observed 95% CL lower limits on the 
contact interaction energy scale A for the electron and muon 
channels, as well as for the combination of those channels. 
Separate results are provided for the different choices of flat 
priors: 1/A 2 and 1/A 4 . 



Channel 


Prior Expected limit (TeV) Observed limit (TeV) 
Constr. Destr. Constr. Destr. 


e + e~ 


1/A 2 
1/A 4 


9.6 
8.9 


9.3 

8.6 


10.1 
9.2 


9.4 
8.6 


/j+fj,- 


1/A 2 
1/A 4 


8.9 
8.3 


8.6 
7.9 


8.0 
7.6 


7.0 
6.7 


Comb. 


1/A 2 
1/A 4 


10.4 
9.6 


10.1 
9.4 


10.2 
9.4 


8.8 
8.4 



CONCLUSIONS 

A search for contact interactions in e + e~ and 
events produced in proton-proton collisions at ^/s = 
7 TeV has been performed. The analysis uses early 2011 
run data amounting to 1.08 (1-21) fb _1 of pp collisions in 
the electron (muon) channel collected with the ATLAS 
detector. The dilepton mass distributions do not display 
significant deviations from the standard model. Using 
a Bayesian approach with a 1/A 2 prior, as was done in 
most previous searches at hadron colliders, the follow- 
ing 95% CL limits are set on the energy scale of contact 
interactions: A" > 10.1 TeV (A+ > 9.4 TeV) in the 
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electron channel and A" > 8.0 TeV (A+ > 7.0 TeV) in 
the muon channel for constructive (destructive) interfer- 
ence in the left-left isoscalar compositeness model. Some- 
what weaker limits are obtained with a prior flat in 1/A 4 . 
These limits are the most stringent to date on /ifiqq con- 
tact interactions and exceed the best existing limits set 
by a single experiment on eeqq contact interactions for 
light-quark flavors. 
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